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fSJ ■ A bijection is given between fixed point free involutions of {1,2, . . . , 2A^} with maximum 

decreasing subsequence size 2p and two classes of vicious (non-intersecting) random walker 
configurations confined to the half line lattice points ^ > 1. In one class of walker con- 

00 ' figurations the maximum displacement of the right most walker is p. Because the scaled 

distribution of the maximum decreasing subsequence size is known to be in the soft edge 
GOE (random real symmetric matrices) universality class, the same holds true for the scaled 

\mJ , distribution of the maximum displacement of the right most walker. 

(-H , Random permutations arc fundamental combinatorial objects, which are intimately related to other 

fundamental combinatorial objects such as Young tableaux via the Robinson-Schensted-Knuth corre- 
spondence. We recall that a Young tableau can be regarded as a numbered diagram of a partition 
Ai > A2 > ■ • • > Ap > 0. The diagram consists of squares drawn within a matrix array with a square 
drawn in each row (1 < j < p) and column k (1 < k < Xj), while in each square is recorded a number 

^ I specified by some rule. Recently random permutations. Young tableaux and their generalizations have 

been shown to be at the core of certain statistical mechanical models of growth processes [E6, 15, 1^, 113, 
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vicious walker paths ||l^, |o[ ||, |T^ and exclusion processes (the latter via mappings to certain growth 
processes and vicious walker paths) amongst other topics. This has led to progress in the study of these 
statistical mechanical models, by way of the progress in the determination of fluctuation formulas for 
^^ I quantities associated with random permutations y, Q, ||] . 

T^ ' As an example of the insight gained, we draw attention to the work of Priihoffer and Spohn |1^ . These 

C^ ' authors identify distinct scaling forms for growth models in the Kardar-Parisi-Zhang (KPZ) universality 

^ ■ class, that is growth models described by the KPZ equation 

C^ ' where ^(i) is a noise term. If the growth profile is curved, the fluctuations are conjectured to coincide 

with the distribution of the largest eigenvalue in the scaled GUE (random Hermitian matrices) , while if 
the growth profile has zero curvature the fluctuations are conjectured to coincide with the distribution 
of the largest eigenvalue in the scaled GOE (random real symmeric matrices). A matrix X from the 
GOE {p = 1) or GUE (/3 = 2) is specified by elements chosen with a joint distribution proportional to 
the Gaussian exp{—/3X^/2). The largest eigenvalue occurs in the neighbourhood of A = ^/2N (which 
is referred to as the soft edge), and by making the scaling A 1-^ \/2N + X/\/2N^^^ the corresponding 
correlation functions have well defined limits |£|. Moreover, in the scaled N -^ 00 limit of both the 
GOE and GUE the distribution of the largest eigenvalue is known exactly in terms of a certain Painleve 
II transcendent pSl p3. This identification was formulated from an exact correspondence between a 
particular model of KPZ growth — the polynuclear growth model — and increasing subsequences of 
random permutations. The latter have been proved to have GUE soft edge fluctuations ||l^, RS, o, nSJ 
in the absence of further constraints, but GOE soft edge fluctuations in the presence of the symmetry 





Figure 1: In diagram (A) there are p = 2 walkers which move a total of 2N = 6 steps. In 
diagram (B) there are N* = 4 walkers which move a total of 2N = 6 steps with the maximum 
displacement of the rightmost walker given hy p = 2. 



constraint restricting the permutations to fixed point free involutions. The two cases correspond to a 
curved and zero curvature interface respectively in the corresponding polynuclear growth model. 

In this work we will identify a statistical mechanical model for which the profile displacement can be 
put into correspondence with the maximum decreasing subsequence length of fixed point free involutions. 
As this quantity has been shown rigorously to have GOE soft edge fluctuations |^, it follows that the 
profile of the statistical mechanical model is in the GOE soft edge universality class. We remark that the 
polynuclear growth model from a flat substrate is also in correspondence with the maximum decreasing 
subsequence length of fixed point free involutions |19[ , the mapping being quite direct (unlike the present 
case). The model to be considered is the random turns model of vicious random walkers. This model can 
be viewed either as a two-dimensional lattice model of non-intersecting directed paths, or as a stochastic 
model of hard core particles on a lattice in one-dimension. In the latter picture, at discrete time intervals 
t = 1,2,... a particle which has a vacant site as its left neighbour or its right neighbour (or both 
neighbours) is selected at random and moved to the vacant neighbouring site (if both sites are vacant, 
either is chosen with equal probability). Plotting the trajectories of the particles on an I — t diagram {I 
labelling the lattice sites) gives the directed, non-intersecting paths picture of the model. 

Our interest is in two classes of configurations of this walker model. The first is when there are 
exactly p walkers, initially equally spaced on neighbouring lattice sites I — 1,2, ... ,p and furthermore 
constrained to the region ^ > 1 (in the vicious walker vernacular, at the site I = there is a cliff at which 
the walkers fall to their death M). It is required that after 2N steps the walkers return to their initial 
sites. 

In the second class of configurations the walkers again begin on the neighbouring sites Z = 1, 2, . . ., 
are confined to the region I > I, and return to their initial sites after 27V steps. But rather than there 
being p walkers there are now N* walkers with N* > N. The parameter p enters by the requirement 
that the right-most walker has a maximum displacement of no more than p lattice sites from its initial 
position. On the other hand the value of N* is not a relevant parameter because only a maximum of N 
consecutive walkers, counted from the right-most walker, move from their initial sites. An example of 
the first and second class of configurations is given in Figure H^. 

Our first result is that both these classes of configurations are in one-to-one correspondence with a 
certain subclass of fixed point free involutions. 



Proposition 1 For both of the two classes of vicious walker configurations specified above there is a 
bijection with fixed point free involutions of {1,2, . . . ,2N} (i.e. permutations consisting soley of two 
cycles) constrained so that the length of the maximum decreasing subsequence is less than or equal to 2p. 

We will give the details of the bijection for the second class of configurations, and afterwards indicate 
the modification required to establish the bijection for the first class of configurations. Denote the walker 
initially at lattice site j by N* + 1 — j {j = 1,2, . . . , N*). A walker configuration can then be coded as a 
sequence of integers from the alphabets 1,2,..., N* and 1,2,..., N* , with the occurence of the integer j 
(j) at position t denoting that walker j (j) moved one site to the right (left) at time step t. 

Of course not all words of length 2N from this alphabet give rise to legal walker configurations. For 
a legal configuration, at each time step we must have that 

ni > n2 > ■ ■ ■ > npf* > (1) 

where 

nj := #j's ~ ^fs 

and after time step 2N, each Uj must equal zero. The requirement (yj) can be represented diagramatically 
as the diagram of a (conjugate) partition in which column j is of length Uj (see Figure 0). A successive 
sequence of diagrams so generated (starting and finishing with the empty diagram 0) uniquely specifies 
the walk thus demonstrating a bijection between such diagrams and walk configurations. 

In the theory of Young tableaux, the diagrams so generated are examples of oscillating tableaux. 
In general these tableaux map onto certain two line arrays p2[ , which in the present case represent 
involutions of {1, 2, . . . , 2A^} with no fixed points. To construct the array, we number the box i if it is 
added to the diagram at step i. If instead a box is removed at step i (say from column j) this is to be 
done via the procedure of reverse column insertion, which means if the particular box ejected, Xi say, was 
then inserted by the Schensted column insertion procedure the original diagram would be restored (see 
[ p^ for a description of the Schensted algorithm) . The fact that the removal occured at step i is recorded 
by putting the pair (i, Xi) into a two line array with i on top. Note that since Xi was bumped out at step 
i, it must have been inserted in an earlier step, so Xi < i. Furthermore all numbers in the array will be 
distinct, and at the end of the procedure there will be N pairs from the numbers {1,2,..., 2A^} with the 
top numbers ordered ii < 12 < • • • < iN- This procedure is illustrated in Figure ^ The pairs forming the 
array can be considered as the two cycles in a fixed point free involution of {1, 2, ... , 2A^}. 

The constraint that the rightmost walker have maximum displacement of exactly p lattice sites to 
the right of its starting position means that the maximum length of the first column of each tableau is 
less than or equal to p boxes. Because each tableaux has the number of the boxes strictly increasing 
down each column and across each row, it follows from the reverse column bumping procedure used to 
form the corresponding two line array that the maximum decreasing subsequence length in the bottom 
line is precisely the maximum size of the first column (see the example of Figure ) which is less than 
or equal to p. Hence the walker configurations are in one-to-one correspondence with the two line arrays 
already noted subject to the additional constraint that the maximum decreasing subsequence length in 
the bottom line is less than or equal to p. In the correspondence between the two line array and fixed point 
free involutions, this constrains the fixed point free involutions to have maximum decreasing subsequence 
length less than or equal to 2p. To see this we note that the fixed point free involution can be constructed 
by extending the top line of the two line array to include all integers 1,2,..., 2A^ in order and filling in 
the bottom line according to pairings implied by the original two line array. We see that if 
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Figure 2: The word corresponding to the walker configuration (B) of Figure ||, the sequence 
of oscillating tableaux corresponding to the word, and the two line array constructed from the 
oscillating tableaux. 



is a particular decreasing subsequence of maximum length q {q < p) in the bottom line of the original 
two line array, then the increasing subsequence of length 2q formed from 



{jl,J2,---,jq}^{Xj,,Xj 



.J 



in the top line of the new two line array gives a decreasing subsequence of length 2q in the bottom line 
of the new two line array. This construction worked in reverse shows that no decreasing subsequence in 
the fixed point free involution can have length greater than 2q. 

The above procedure associating each walker configuration with a two line array is reversible in that 
starting with a two line array of the type specified a unique sequence of oscillating tableaux and thus 
walker configuration can be constructed. Following [ p2| we work backwards in the construction of the 
two line array from the sequence of oscillating tableaux. In going from the tableau at step i to that at 
step i — 1 there are two distinct situations. One is that i does not appear in the top row of the two line 
array, indicating that the tableau at step i was not the result of removing a box from the tableau at step 
i — 1, but rather came from adding a box labelled i to the tableau at step i — 1. Thus deleting the box 
labelled i from the tableau at step i gives the tableau at step i — 1. On the other hand we may have 
that i does appear in the top row of the array, being part of the pair {i,Xi). In this case the tableau 
at step i was obtained from the tableau at step i — 1 as a result of an inverse column bumping which 
ejected Xi. Thus the tableau at step z — 1 is constructed from the tableau at step i by Schensted column 
inserting Xi. An example of this inverse procedure is given in Figure y[ From the rules of the column 
insertion the maximum attained length of the first column of the tableaux will equal the length of the 
largest decreasing subsequence in the bottom line of the two-line array and thus be less than or equal to 

P- 

Hence for every sequence of oscillating tableaux, starting and finishing with the empty tableau and 
having column length less than or equal to p, there is a two line array equivalent to a fixed point free 
involution having maximum decreasing subsequence length less than or equal to 2p, and furthermore 
the correspondence can be established in the reverse direction. Because there is a bijection between the 
oscillating tableaux and random walker configurations, the result of the Proposition |l| for the second class 
of walker configurations is established. 

Let us now turn our attention to the first class of configurations. The walker configurations are again 
written as words, this time from the alphabets 1,2, ... ,p and 1,2, ... ,p. The constraint (^ (with N* 
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Figure 3: Correspondence between a two line array corresponding to a fixed point free involution 
and a sequence of oscillating tableaux, and the correspondence between the oscillating tableaux 
and a word. The word is equivalent to walker configuration (A) of Figure |^, but translated at 
least one lattice site to the right so that N* > A^, with stationary walkers filling the intervening 
sites to the left down to / = 1. 
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Figure 4: The oscillating tableaux for the configuration (A) of Figure |^, or equivalently the 
word of Figure 0, constructed from the rules for a fixed p of walkers, and the corresponding two 
line array. 



replaced by p) is represented as a diagram but now with row j of length Uj rather than column j as 
previously, this feature being the essential difference between the two cases. Note that the length of the 
first column now represents the number of walkers displaced from their initial conditions. 

As before, the boxes are numbered by i if added at time step i, and removed via the reverse column 
bumping procedure, with the fact that the removal occured at step i recorded by putting the pair {i,Xi) 
into a two line array with i on top. An example is given in Figure 0. The constraint that there be less 
than or equal to p walkers restricts the first column length to be less than or equal to p. As already 
noted, the fact that the reverse column bumping procedure is used to construct the two line array from 
the tableau implies the former must therefore have decreasing subsequence length no greater than p. 
Thus each walker configuration can be mapped to a unique two line array of the same type as occured 
in the corresponding mapping for the second class of configurations. Futhermore, we have detailed how 
to associate such two line arrays with a unique sequence of oscillating tableaux. From this sequence of 
oscillating tableaux we can construct the word corresponding to the walker configuration. Note that this 
differs from the construction in the case of the second class of configurations because now it is row j which 
specifies the moves of walker j. The final result is that there is a bijection between the first class of walker 
configurations and two line arrays with top line ordered «i < 12 < • • • < iv: ^^rid maximum decreasing 
subsequence length no greater than p. The latter being equivalent to fixed point free involutions of 
{1,2,..., 2A^} with maximum decreasing subsequence length no greater than 2p, we see that Proposition 
|l| is now established. 



The second class of configurations count the number of walker configurations with a specific bound p 
on the maximum displacement of the right- most walker. From the / — t diagrams of Figure n^ we see that 
this is equivalent to counting the number of configurations which give rise to a growth profile with a bound 
p on its maximum spread, Ln say. Our interest is in the distribution of L^. Now, with Lt denoting the 
displacement of the right- most walker after t steps, the symmetry of the configurations under t -^ 2N — t 
means the statistical properties of Lt are the same as those of L2N-t- At the centre of symmetry will be 
the maximum displacement Ln, and the conjecture of Prahoffer and Spohn predicts GUE fiuctuations if 
the profile is curved at this point, or GOE fiuctuations if the profile has zero curvature. Unfortunately 
the analytic form of the profile is not known, so we cannot make use of this prediction presently. 

In fact the nature of the fluctuations can be determined rigorously by using the bijection of Proposition 
H between the second class of configurations and fixed point free involutions. In particular we can make use 
of the known distribution of the maximum decreasing subsequence length for fixed point free involutions 
of {1,2,..., 2N} to deduce the limiting distribution of L^. Regarding the former, let U§^ denote the 
maximum decreasing subsequence length, and define the scaled quantity 



.inv , 


L'^^ - 2V2N U^^/2 ~ 2^N/2 


.N ■ 


(V27V)i/6 (^iV/2)i/6 




lim Prfxr < a;) = Fi{x) 



Then it is proved in [pj that 

where Fi{x) denotes the cumulative distribution of the largest eigenvalue of matrices from the scaled 
GOE |2J, |l^ . The following result is then an immediate consequence of Proposition |l|. 

Proposition 2 Let Ljv denote the maximum displacem.ent of the right-most walker in the second class 
of random walker configurations specified above, and set 



Ln - V2N 

XN 



i(27V)i/6 ■ 

Then 

lim Pr(xjv < a;) = Fi{x). 

Hence the walker profile at its maximum width exhibits GOE fluctuations. The converse of the prediction 
of Prahoffer and Spohn would then imply that the walker profile has zero-curvature at this point. 

As a final issue we consider the p-dimensional integral formula for the number, fj^^ say, of fixed point 
free involutions of {1, 2, ... , 2A'^} constrained so that the length of the maximum decreasing subsequence 
is less than or equal to 2p. With USp{p) denoting the group of 2p x 2p unitary symplectic matrices (or 
equivalently the group oi p x p unitary matrices with real quaternion elements ) , it was shown by Rains 
[Ml that 



^ ) P- -^0 Jo . 1 _;^ l<J<fe<P 



where Zj = e^ ^ . 

The formula (0) is in fact a special case of a counting formula for a class of vicious walker paths. Thus 
consider p vicious walkers in the lock step model, confined to the lattice sites I > I, starting at positions 

1<4°^< 4°^ <...<;("), (3) 
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and arriving at positions 

1 < ^1 < ^2 • • • < ^p (4) 

after n steps. With Zn{li , ■ ■ ■ , Ip , h, ■ ■ ■ , Ip) denoting the number of distinct walker configurations of 
this prescription, we have the foUowing result. 

Proposition 3 

7 (l'^°'> /(O)- 7 7 "1 

■^n V"! J ■ ■ ■ J ""p J '■li • ■ ■ I 'pj 

ddi- 
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d9p(^2 cos Bj ] det 



Ah-i'-°^)ej _ p^ih+i'-°^)0, 



j=i 



j,k=l,...,p 



(5) 



Analogous to the proof of a similar counting formula in B , this can be verified by first noting from 
the definition of the particular lock-step model that 

■^n(,^l : ■ • ■ T ^p 'i ''1: • ■ • : ''p) -^ -^nv^li ■ ■ • : ''p) 

is the unique solution of the multidimensional difference equation 

■Z^ra+l(^l, ■ • ■ )^p) = Zn{ll ~ Ijh, ■ ■ ■ ,lp) + Znihjh — ^, ■ ■ ■ ,lp) 
+ • • • + Zn{ll, h, ■ ■ ■ ,lp — ^) 
+ Zn{ll + 1, Z2, • ■ • , ip) + Zn{ll, Z2 + 1, • . • ,^p) 

+ --- + z„(/i,;2,...,;p + i) (6) 

subject to the non-intersection condition 

z„Gi,...,;p) = o if ij^ikij^k) (7) 

the constraint Ij > 1 (j — 1, . . . ,p) which requires 

Znill,...,lp)=-0 if h^O (8) 

(here use has been make of the ordering (||)) and the initial condition 



z„(z(°\...,;(o);/i,...,g = n<5, 



I ",h 
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where again use has been make of the orderings (0) and 
To verify that @ satisfies (||) we note that (||) gives 
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J2=Ai + l P 
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and recalling (|5|) we can immediately identify the right hand side of (|9|) with the right hand side of (^. 
To verify (R) we simply note that if Ij = Ik for any j ^ k then two rows of the matrix in dq) are the 
same so the determinant vanishes. The condition (W is a property of (0) since with ^i = the integrand 
is odd in 9i and thus the integral vanishes. Finally, to verify the initial condition (Qj we make use of the 

definition of a determinant 

p 

det[ajk]j,k=i,...,p ^ ^ e{P)Y[ajP{j), 
PeSp ]=i 

where e{P) denotes the parity of the permutation P, to expand the integrand in dq) and integrate term 
by term. Recalling each I, and h is positive, this gives 



Zo{l^°\ . . . , ;(0); Zi, . . . , Zp) = ^ e{P) f[ <5,^^,(o,^. 



(11) 
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The ordering constraints (0) and (0) imply that all terms in (0) except for the identity permutation 
must vanish, and so (^ is indeed satisfied. 

Since the difference equation, the boundary conditions and the initial conditions are all satisfied by 
(||), we conclude that (||) correctly represents Zn{h, ■ ■ ■ , Ip)- 

Although Zn is a positive integer, the integrand in dq) is complex. An integral representation with a 
positive real integrand can, in the case l-=lj{j — l,...,p), be obtained by first noting 



^ni+n2 V'l I ■ ■ ■ J 'p I 'l I ■ ■ ■ J 'p ) 

/ J -^nii'l I ■ ■ ■ I 'p I '1: ■ ■ ■ J fpj-'irig I'l: ■ ■ ■ I 'pi 'l : ■ ■ ■ i 'p )■ 

l<h<l2<-<lp 

From (p[) , and after simple manipulation of the determinant therein, we see that 

Zni\li , • • ■ , tp ; '1, • • ■ , 'pj'^n2 v'l? • ■ ■ I 'pi '1 7 • ■ ■ J 'p ) 

—^ dOi--- dOp d4>i--- d(^p(^2cos0jj (^2cos(^ 
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j,k=l,...,p 






e "*= 



^^i«)0, 



i,fe=i, 



(13) 



This is an even symmetric function of the Z^'s, which vanish for Ij — Ik and Ij — 0. Consequently the 
sum in (03) can be replaced by 

-. oc oc oc 



2Pp\ 



11— — 0012 — — OQ Ip — — 00 



Performing the sum in ( |l3| ) using 

00 

/ — — C30 

and putting ni + n2 = n gives 

Z (';(o) /(o)./(o) /(o)N 

■^nA'l : ■ ■ ■ I 'p I 'l I ■ ■ ■ I 'p y 

{27T)P2Pp\j_, 7_^ n^ 



2 cos 00 



det 



e fc J — e * 



j,k=l,...,p 



(14) 



We can use (^4|) to rederive (||) since from the definitions we have 

fNp^Z2Nili ,.. .,/^°Vl,- •■ ,^p) ,(0) , .,. , ,• (15) 



if>=h=jU=i,...,p) 



;(0) _ 



Setting l-'=j{j = l,...,p)m (|14|), noting that with Zj = e ^ we have from the type C Vandermonde 
formula |2^ 

det[z^^-z7%-,,.i,...,,p = n|l-zfP n \i~^,^k\'\zj-z,\', (16) 

j=l l<j<k<p 

and making use of the fact that the integrand is even we see that the formula (0) indeed results. 

Acknowledgements 

PJF acknowledges the financial support of the Australian Research Council, and that of funds obtained 
by Prof. K.Aomoto for his visit to Japan in June 2000 which facilitated the present collaboration. 

References 

[1] J. Baik, P. Dieft, and K. Johansson. On the distribution of the length of the longest increasing 
subsequence of random permutations. J. Amer. Math. Soc, 12:1119-1178,1999. 



[2] 
[3] 

[4] 
[5] 
[6] 

[7] 
[9 

[lo; 

[11 

[12 
[13 



J. Baik and E.M. Rains. Algebraic aspects of increasing subsequences, math. CO/990508"^, 1999 



J. Baik and E.M. Rains. The asymptotics of monotone subsequences of involutions. 



math.CO/9905084, 1999 



J. Baik and E.M. Rains. Symmetrized random permutations. math.CO/991001£ , 1999 



J. Baik. Random vicious walks and random matrices. Comm. Pure Appl. Math., 53:1385-1410, 2000. 

A. Borodin, A. Okounkov, and G. Olshanski. Asymptotics of Plancherel measures for symmetric 
groups. J. Am. Math. Soc, 13:481-515, 1999. 

P.J. Forrester. Probability of survival for vicious walkers near a cliff. J. Phys. A, 22:L609-L613, 
1989. 

P.J. Forrester. Finite size corrections to the free energy of Coulomb systems with a periodic boundary 
condition. J. Stat. Phys., 63:491-504, 1991. 

P.J. Forrester. The spectrum edge of random matrix ensembles. Nucl. Phys. B, 402:709-728, 1993. 



P.J. Forrester. Random walks and random permutations, math. CO/9907037 , 1999 



P.J. Forrester. Painleve transcendent evaluation of the scaled distribution of the smallest eigenvalue 



in the Laguerre orthogonal and symplectic ensembles, alin. SI/0005064 , 2000 



W. Fulton. Young Tableaux. London Mathematical Society Student Texts. CUP, Cambridge, 1997. 

A.J. Guttmann, A.L. Owczarek, and X.G. Viennot. Vicious walkers and Young tableaux i: without 
walls. J. Phys. A, 31:8123-8135, 1998. 



[14] C.A. Tracy J. Gravner and H. Widom. Limit theorems for height fluctuations in a class of discrete 
space and time fluctuations, [math. PR/0005133, 2000. 



[15] K. Johansson. Discrete orthogonal polynomial ensembles and the Plancherel measure. 



math.CO/9906120, 1999 



[16; 

[17; 
[is; 

[19 

[20 

[21 

[22; 
[23; 



K. Johansson. Shape fluctuations and random matrices. Commun. Math. Phys., 209:437-476, 2000. 

C. Krattenthaler, A.J. Guttmann, and X.G. Viennot. Vicious walkers, friendly walkers and Young 
tableaux: 11. J. Phys. A, 33:8835-8866, 2000. 

A. Okounkov. Random matrices and random permutations. Int. Math. Res. Notices, (20): 1043-1095, 
(2000). 

M. Prahofer and H. Spohn. Universal distributions for growth processes in 1 + 1 dimensions and 
random matrices. Phys. Rev. Lett, 84:4882-4885, 2000. 

R.A. Proctor. Odd symplcctic groups. Inv. Math., 92:307-332, 1988. 

E.M. Rains. Increasing subsequences and the classical groups. Elect. J. of Combinatorics, 5:#R12, 
1998. 

S. Sundaram. The Cauchy identity for Sp{2n). J. Comhm. Theory A, 53:209-238, 1990. 

C.A. Tracy and H. Widom. Level-spacing distributions and the Airy kernel. Commun. Math. Phys., 
159:151-174, 1994. 



[24] C.A. Tracy and H. Widom. On orthogonal and symplectic matrix ensembles. Commun. Math. Phys., 
177:727-754, 1996. 



10 



